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Abstract 

When all Green's functions are known in a given gauge we may raise a question of 
whether it is possible or not to derive the corresponding ones in a different gauge. The 
answer is negative in QCD but affirmative in QED provided that we confine ourselves 
to the covariant gauge characterized by a gauge parameter a. We shall discuss the 
physical significance of this conclusion. 



1 Introduction 

In quantizing a gauge-invariant Lagrangian we encounter a well-known difficulty in finding 
the canonical conjugate of the time-component of the gauge field. This is a reflection of the 
non-uniqueness of the solution of the gauge field equation due to the gauge freedom. This 
difficulty has been resolved in QED by Fermi's introduction of the gauge-fixing term into 
the Lagrangian that provides us with the lacking canonical conjugate. 

Once the theory is quantized by this method, however, the introduction of indefinite 
metric is indispensable since the vector field is obliged to inherit it from the Minkowski metric. 
This means that the state vector space resulting from the quantization of the gauge field is 
larger than needed for physical interpretation, and we have to pick out its physical subspace 
by introducing a subsidiary condition. This condition plays a dual role of eliminating the 
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indefinite metric inherent in the state vector space as well as of recovering the classical gauge 
field equation, not modified by the gauge-fixing term, in the physical subspace. 

In the present paper we shall confine ourselves to the so-called covariant gauge specified by 
a parameter a called the gauge parameter. Quantized gauge theories are no longer invariant 
under local gauge transformations because of the presence of the gauge-fixing term, but it so 
happens that the resulting theory is invariant under new global transformations called the 
BRS transformations [1]. 

In Section 2 we shall reinstate the essence of the BRS transformations in connection 
wit the gauge-dependence of Green's functions and shall clarify the condition under which 
Green's functions become gauge-independent. In tackling the problem of gauge-dependence 
the renormalization group (KG) approach [21 [3] is useful and it is briefly recapitulated in 
Section 3 with an emphasis on the gauge field propagator. 

The question of gauge-dependence was once raised by Shirkov in connection with the 
use of RG in perturbative QCD [4J. The results depend sensitively on the renormalization 
scheme so that we propose here a scheme of making the beta-function independent of the 
gauge parameter a. 

In Section 4 the RG equations for the running coupling constant and gauge parameter 
are solved and their ultra-violet asymptotic limits are studied to answer the question of 
gauge-dependence of Green's functions. For this purpose we derive a sum rule [5l [6] which 
enables us to express the renormalization constant Z^^ in a simple form. 

Green's functions defined as the vacuum expectation values of the time-ordered products 
of BRS invariant operators are always independent of a, but those made up of BRS variant 
ones are a-dependent. For the latter an important question is whether it is possible or not 
to continue Green's functions analytically as functions of a. When it is possible to continue 
them from a\ to a2, we say that a\ and 02 are connected. Then, mutually connected values 
of a form a set called an equivalence class of gauges, and we may ask how many equivalence 
classes of gauges there are in QCD and also in QED. The answer is that three in QCD and 
one in QED, respectively. QED is simple, but QCD is complicated and we may ask what the 
multiplicity of equivalence classes would mean. We shall discuss the physical significance of 
these results in connection with the gluon mass. 

In Section 5 this quest is further pursued by extending the gauge parameter a into the 
complex plane. 
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2 BRS Invariance 



Local gauge transformations in classical gauge theory are replaced by global BRS transfor- 
mations and we shall briefly recapitulate their properties 



BRS transformations 

The standard Lagrangian density of a gauge theory, say QCD, is given by 

>C = Cinv + J^gf + J^FP (2-1) 

where Cinv denotes the classical gauge- invariant part, Cgf the gauge-fixing terms and jOfp 
the Faddeev-Popov (FP) ghost term characteristic of non-Abelian gauge theories: 

J^inv = —Ff,^ ■ F^u -'4'ilnDi^ + m)ilj , 
Cgf = d^B-A^ + ^aB-B, 

Cfp = id^c-D^c. (2.2) 

in the customary notation. The gauge parameter is denoted by a and represents the 
covariant derivative whose explicit forms are given by 

D^ij = {d^-igT-A^)i;, 

D^c = d^c + gA^xc. (2.3) 

The BRS transformations of the gauge field A^ and the quark field t/j are defined by 
replacing the infinitesimal gauge function by the FP ghost field c or c in their respective 
infinitesimal gauge transformations: 

SAi^ = Df,c, Sijj = ig{c-T)'ijj , 

= D^c, SiP = zg{c-T)'ip. (2.4) 

For the auxiliary fields B, c and c we require 

SC^SC^O, (2.5) 

then we find 

5B = 0, 5c = iB, 5c=~g{cxc), 

SB^O, Sc^iB, Sc^~g{cxc), (2.6) 

where B is defined by 

S + S -i^(c X c) = 0. (2.7) 
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In general the BRS transforms of a field are given in terms of the BRS charges Q b and 
Qb by 

5cj) = t[QBA]T^ 5^ = i[QbA]t^ (2-8) 
QI = QI = QbQb + QbQb = ^. (2.9) 

We choose the — (+) sign in (12.91) when is even (odd) in the ghost fields c and c that are 
anticommuting hermitian scalar fields. 

The sum of the gauge-fixing and the FP ghost terms can be expressed as 

Cgf + CFP = 5{-id^c-A^-\c-B) (2.10) 

and evidently we have 

5A™ = 0. (2.11) 
Namely, Cinv is closed and Cgf + Cpp is exact, and 

5£ = 0. (2.12) 



BRS Cohomology ^\ 

The quantization of the gauge field and the introduction of the auxiliary fields B, c and 
c introduce indefinite metric into the state vector space V. 

A physical state |/) is defined by the following subsidiary condition: 

Qb|/) = 0, |/)gV. (2.13) 

In particular, the vacuum state |0) is physical, 

Qb\0) = 0. (2.14) 

The physical subspace Vphys is then defined by 

VpHys = {\f)\ Qb|/) = 0, |/)GV}. (2.15) 

It is essentially a collection of closed states with respect to the nilpotent operator Qb- We 
also introduce a subspace Vd defined by 

Vd = {QB\f)\ |/)GV}. (2.16) 

This is a collection of exact states with respect to and the Hilbert space TC is defined as 
the BRS cohomology by 

n = Vphys /Vd. (2.17) 
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Then let us consider a set of closed operators A, B, - ■ ■ satisfying 

SA = 6B = --- = 0, (2.18) 

then 

{0\6M -AB-- -10) = {0\6{MAB ■ ■ ■ |0) = , (2.19) 

since the vacuum state is physical. 

Let Cj and £// be two BRS invariant Lagrangian densities, namely, 

SCi = SCn = . (2.20) 

Furthermore, let us assume that their difference be exact: 

AC = 6Cii - 6Ci = 6M . (2.21) 

For instance, Lagrangian densities corresponding to two distinct values of a in (12.21) satisfy 
these two conditions since we find 

M=-hAa{c-B). (2.22) 

Then we introduce time-ordered Green's functions in two gauges given above. 
In the gauge I we have the following path integrals: 

{AB---)i = ^ I V{path) AB--- exp{iSi) (2.23) 
Ni J 



Ni = j V{path) exp{iSi) . (2.24) 



V{path) denotes the path integral over all the field variables, and we have similar expressions 
in the gauge II, and the difference between the two actions is given by 

AS = Sn-Si = j d^x ^'^ = j . (2.25) 

Green's functions in the gauge II can be expressed as 

/AR \ {AB ---expjt AS)) J 

{AB---}jj = . 2.26 

{exp{tAS)}j 

Now we expand the denominator on the r.h.s. of (12.261) in powers of AS* and use (I2.2ip and 
(EUD to obtain 

{exp{iAS))i = 1. (2.27) 

Thus (12.261) reduces to 

{AB ---)n = {AB - - - exp{i AS))i . (2.28) 
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The r.h.s. can be expanded in powers of Aa = aji — aj. When this expansion converges 
this is an analytic continuation of Green's functions as functions of a. When this is the case 
we may say that the two values of a, ai and 02, are connected and the set of a mutually 
connected will be called an equivalence class of gauges. The question of how many classes 
there are in a given gauge theory will be discussed in a later section. 

When all the operators are closed, satisfying (12.181) . we obtain, with the help of (12.191) 
the equality 

{AB---)n={AB ■■■)!. (2.29) 

This shows that Green's functions cosntructed in terms of closed operators alone are gauge- 
independent. The S'-matrix elements for observable hadronic processes are obtained by 
applying the LSZ reduction formula [7] to Green's functions defined in terms of BRS invariant 
operators so that they are independent of the choice of the gauge parameter a. Of course, 
this statement does not apply to Green's functions constructed in terms of BRS variant 
operators, such as the gauge field propagator which we shall investigate in a later section. 

3 Renormalization Group 

Once Shirkov emphasized the gauge. dependence of the KG treatment in perturbative QCD [4j 
since the results are sensitive to the renormalization scheme. He even presented an example 
in which asymptotic freedom is valid only for positive values of a but not for negative ones. 
This is apparently due to the a dependence of the beta-function. Therefore, we shall first 
briefly review a renormalization scheme which makes the beta-function independent of the 
gauge parameter [6]. 

First we shall refer to ( 12.28^ which has been derived in the unrenormalized version. The 
field operators Ax and B and the gauge parameter a are multiplicatively renormalized: 

= Zl^\a) Ax , = Z-^^^{a) B , «(°) = Z^ia) a , (3.1) 

where the superscript (0) is attached to the unrenormalized expressions. 

As a special case of (12.281) we shall choose an = a and «/ = 0, then we have 

{AB---)a = (^AB---exp(^-^ j d^xB{x) ■ B{x)^^ , (3.2) 

in the unrenormalized version. The renormalized version of (13.21) takes exactly the same 
form because of the identity 

c,{o)5(o).5{o)^^5.5_ ^33) 
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Then we introduce RG equations in the Landau gauge, a = 0. For the renormahzed version 
of Green's functions we have an equation of the following form: 

(I)o + 7g)(^5---)o = 0, (3.4) 

where 7g is the anomalous dimension of Green's function {AB ■ • ■ )o in the Landau gauge 
and 

^^o = /i|^ + /3(^7)^. (3.5) 

The anomalous dimension of the gauge field in the Landau gauge, which is obviously a- 
independent, is denoted by ■jv 

Then the l.h.s. of (13. 2p . renormalized in the Landau gauge, satisfies 

iVo + ^G){AB---)^ = 2^v(AB---exp('-^ [ d'^x B{x) ■ Bix)]) . (3.6) 
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On the other hand, we also have 



d 



la 



a—{AB---)^ = (^AB---exp\^—jd^xB{x)-B{x)jJ . (3.7) 

Combining these two equations we find 

(2) + 7g)(A5---)„ = 0, (3.8) 

where 

V = Vo + 2a7y^ • (3.9) 
Next we introduce the gauge field propagator 

(^a(^),4(?/)) = (^^'^^ / d'ke^'^^^-y^D,,{k) , (3.10) 

which is the vacuum expectation value of the time-ordered product of two color gauge field 
operators, and 

- (^--^^)^(*^") + "^. (3,11) 

D(k\a) = ! dm \, f'"'' . . (3.12) 



kP' + Tm? — ie 
Then introduce 

R{k'^,a) = P D{P,a) . (3.13) 
Since the gauge field has been renormalized in the Landau gauge, we have 

i?(/i^ 0) = 1. (3.14) 
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In other gauges we have to employ the renormahzation factor ^3(0) instead of ^3(0), and 
this amounts to a further renormahzation for R. The properly renormalized R function and 
the gauge parameter are denoted by R and a, respectively, and are given, in a consistent 
manner with (13. ip . by 

R{k^,a) = R{fj,'^,a)R{k'^,a) , (3.15) 
a = R{i/,a)a. (3.16) 

The function R{k^, a) renormalized in the Landau gauge satisfies 

{V + 2-fv)R{k\a) = 0. (3.17) 

Then the function R{k^, a) normalized by 

R{f/,a) = l (3.18) 

satisfies 

(D + 27y)^(A;2, a) = , (3.19) 
where the new anomalous dimension 7^ is given by 

^v = lv + lv\nR{fj.\a). (3.20) 

Next we switch the set of parameters from (/i, g, a) to (/x, g, a). The anomalous dimension 
7y will be expressed as a function of g and a, and it will be denoted by •jviOy C() from now 
on. Then, 

= /i^ + /?(^7)^-2a7y(^,«)— . (3.21) 

In this way we have established a renormahzation prescription leading to an a-independent 
beta-function. For the quark field we can derive its anomalous dimension 7^ in a manner 
similar to the above derivation of jv Since the beta-function is independent of a, the 
concept of asymptotic freedom is gauge-independent. From now on we shall skip the bars 
introduced above. 



4 Asymptotic Limits of Running Parameters 

The RG equations for QCD have been studied in detail and we shall reinstate their essence 
in what follows ^5, 6J. 
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Renormalization Constants 

An element of RG may be expressed as 

R{p) = expipV) , (4.1 

where p denotes the parameter of RG and the composition law of this group is given by 

Rip)Rip') = Rip + p') . (4.2 

Let Q be a function of g, a and p, and we define the running Q to be 

Q{p) = exp{pV)Q{g, a, p) = Q{g{p), a{p),p{p)) , (4.3 

with the initial condition 

Q{0) = Q. (4.4 
Let the anomalous dimension of Green's function G{pi; g, a, p) be 7(5', a), then we have 



[1^ + lig, Oi)]G{pi; g,a,p) = 0. 
Its running version is defined by 

G{p) = exp{pV) ■ G{pi] g, a, p) , 

and it satisfies 



d_ 

dp 

The formal solution is given by 

G{pi;g,a,p) = exp 



G{p) = -^{p)G{p) . 



dpW 



■ G{pi]g{p),a{p),p{p)) . 



(4.5 

(4.6 
(4.7 

(4.8 



Then we assume, in the presence of a cut-off A, that the running coupling constant g{p 
tends to the unrenormalized or the bare one go in the limit p — 00, namely, 



(0) 



and similarly 



lim g{p) = g 



lim a{p) = a^^^ , 

p— >oo 

lim p{p) = p lim = 00 . 

p—fCX) p-^oo 



In this limit (14.81) reduces to 

G{pi;g,a,p) = exp 



dp' lip 



{Pi;g'"\ a'"', 00} 



(4.9 

(4.10 
(4.11 

(4.12) 
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where G^^^ denotes the unrenormahzed version of Green's function G. Then the renormal- 
ization constant of G denoted by Z is given by 



Z = exp 
The solution of (13. 190 reads as 
R{k'^; g, a, /i) = exp 



dpw: 



dp'lvip) 



R{k ■,g{p),a{p),p{p)) . 



exp 



dp'lvip') 



(4.13) 



(4.14) 



Now apply the Lehmann representation (13.121) to the l.h.s. and take the limit p ^ oo after 
putting k"^ = /i^(p), then with the help of (14.131) we obtain 



(4.15) 



In the cut-off theory we first take the limit p —>■ oo and then A ^ oo, but in what follows 
we invert the order of limiting procedures by taking the limit A — > oo first. Thus some 
of the initial conditions introduced in the cut-off theory are not necessarily satisfied. As 
an example we shall see later that the limiting values (14.91) and (14.101) cannot be arbitrary 
despite our expectation that the unrenormahzed g^^^ and a^^^ should be chosen arbitrarily. 



Asymptotic Limits of Running Parameters 

Next we shall study the RG equations for the running parameters g{p), a{p) and /2(p) 
that follow from (14.31) . 

^5(p) = (4.16) 

Att(p) = -2a{p)ivip), 
j-^p{p) = Hp) or Ji{p)=peP. 
First we shall define their asymptotic limits by 

g{oo) = goo, a(oo) = ctoo, /i(oo) = oo . (4.17) 

In the absence of a cut-off they do not necessarily reduce to their unrenormahzed counterparts 
in dUD and fICTD . 

In perturbation theory (3{g) and ivig, ct) are given in the form of power series, 

Pig) = g%Po + Pig' + ■■■), (4.18) 

7y((7,a) = /(7o(«) + 7i(«)/ + ■■■), (4-19) 
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where 



7o(a) = 700 + 701", 

7i(«) = 7io + 7ii" + 7i2«^ (4.20) 



The lowest order coefficients are given by 



When is negative, namely, when Nj < 16, asymptotic freedom is realized and we shall 
assume it in what follows. 

Asymptotic freedom [H [9] is characterized by 



0, (4.22) 



and for large values of p we obtain approximately 

1 



f(p) = ^ (& = -2/3o>0). (4.23) 



By integrating the second equation in fl4.16p we find a sum rule [H [6]: 



a 



00 



In— = -2/ dpMp), (4.24) 
a 



or 

r /"oo 1 „ 



h = exp 



dplvip) 



a 



(4.25) 



00 



The essence of our argument is based on this sum rule. 

Our main problem is the determination of the asymptotic limit a^o, but this problem has 
been discussed in detail before, so we shall only quote the results in what follows. We shall 
come back to it, however, in the next section for a different purpose. 

For a = we have a{p) = and hence ctoo = 0. It is clear that a{p) and a are always of 
the same signature. We find three possibilities for aoo- 

aoo = -00, 0, ao , (4.26) 

where ao is defined by 

7o(ao) = 7oo + 7oiao = or ao = -700/701 • (4.27) 
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For ao > OT Nf < 10, we find 

aoo = <( 0, a = (4.28) 




and in this case the sum rule can be expressed as 



/ 



dm" p{m^) = —e{a) (4.29) 
where 

= (4.30) 
10, a < . 

The r.h.s. of fl4.29p is a continuous function of a, but its derivative with respect to a 
develops a discontinuity at a = 0. Thus, the set of real a, denoted by D, is divided into 
three connected subsets: 

D = D(-oo) IJd(O) IjD(ao) , (4.31) 
where we confine ourselves to positive ao or Nj < 10 

D = {a\ a G -R} , 

D{-oo) = {a\ a < 0} , 

D(0) = {a\ a = 0}, 

L)(ao) = {a\ a > 0}. (4.32) 

Thus, for gauge-dependent Green's functions such as the gluon propagator we find three 
equivalence classes of gauges. 

In one of the papers by the present authors [lO] the residue of the massless pole in the 
two-point function 

{Al{x),FlM) (4.33) 

has been studied, and it has been concluded that the residue vanishes in two classes D{—oo) 
and -D(O), suggesting that the gluon in these gauges would likely be massive, whereas it is 
non- vanishing in the class -D(ao) implying zero mass for the gluon. This means that the 
gluon mass would be gauge-dependent. The interpretation of this result is delicate. In that 
paper we tacitly assumed that the gluon mass is a physical quantity and that the above 
result indicates the gauge-dependence of QCD, or more precisely the class-dependence of 
the gluon mass, namely, it is zero in the equivalence class of gauges D(ao), but it is non-zero 
in the classes D{—oo) |JD(0). The gluon mass denoted by M must be RG invariant 

VM = (4.34) 
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within a class. 

We found that this interpretation is misleading since the physically observable quantities 
must be gauge-independent as illustrated by (I2.29p . 

The S-matrix elements for hadronic processes are gauge-independent as is clear from 
f l2.29p , and the condition for color confinement [5l El [HI [121 IE] 

Zg-^ = (4.35) 

is satisfied in the classes D{—oo)[J D{0). Therefore, in these classes color confinement is 
realized and the unitarity condition of the S*- matrix between two hadronic states \a) and \b) 
reads as 

J2{b\S^n){n\S\a) = {b\a) , (4.36) 

n 

where the sum over the intermediate states is saturated by hadronic states alone without 
introduction of the confined quarks and gluons. Since the S'-matrix is gauge- independent 
this statement is also valid for the class D{ao). Thus, confinement is a gauge- independent or 
class-independent concept. This conclusion does not contradict the gauge- dependence of the 
gluon mass since the mass of the confined gluons is never observed and turns out to be an 
unphysical quantity. A similar observation has been made by Fujikawa, Lee and Sanda p3] . 
that particles with gauge-dependent masses are unphysical and not subject to observation 
in connection with the gauge. 



5 Introduction of the Complex Gauge Parameter 

From (14.191) and (I4.20p we can readily deduce that 7y can be expanded into a double power 
series in and ag^. For sufficiently large p, we can assume that g'^{p) ^ 1 and if a{p) is 
bounded below a certain constant, we may use perturbation theory in powers of g^. In the 
lowest order we have 

f^ = m-^Pof = ~f, (5.1) 



and its solution is given by 



When a is bounded we employ the following approximate equation: 

do. 

— = -2a^v ~ /(p)a(ao - «) , (5.3) 
dp 

where 

f{p)=2-fo,f{p)>0. (5.4) 
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When the integral (14.241) is convergent, Ooo niust be finite. The expansion of 7y in powers 
of g'^ starts from g'^ and for large values of p the behavior of g'^{p) is given by (14.231) so that 
the integral of the first term in the expansion of 7v diverges like 



oo 





(5.5) 



In order for the integral of the power series to converge term by term in (14.241) . therefore, 
the condition 

7o(«oo) = (5.6) 
must be satisfied. With reference to (I4.27P we find 

"oo = "0 (5.7) 
When the integral (I4.24p is divergent, ln(a/aoo) should diverge so that we find 



Ono = or ± oo . 



(5.^ 



Now, when a goes off aoo starting from its neighborhood for increasing p, is called 
a repulsive asymptotic limit. Otherwise, when a approaches a^o again starting from its 
neighborhood, it is called an attractive asymptotic limit. 



Case of positive ao {Nj < 10) 

Integration of (15.31) in the neighborhood of ao yields 



"0 - 



exp 



-ao / dp' j\p') 
Jo 



, for p ^ oo - 



ao — a 

Thus ao is found to be attractive. In the neighborhood of 0, on the other hand, we find 



(5.9) 



a(p) 



a 



exp 



ao / dp'fip') 
Jo 



oo , for p ^ oo 



(5.10) 



Thus aoo = is found to be repulsive. 

The fiow of a(p) on the real a axis for increasing p is given in Fig. 1. 



-oo 



+ CX3 



Figure 1: Flow of a(p) along the real a axis for increasing p (ao > 0). 



From Fig. 1 we conclude that — oo is also attractive. 
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< \ > \ < 

-oo q/q +00 



Figure 2: Flow of a(p) along the real a axis for increasing p (ao < 0). 

Case of negative ao (10 < Nf < 16) 

In this case and —00 are attractive and ao is repulsive. 

In what follows we shall confine ourselves to the case of positive ao and shall study what 
would take place when a turns out to be complex while keeping real. For this purpose 
we assume that g'^ is already sufficiently small and that |a(p)| is bounded to guarantee the 
following treatment of the RG equations. 

The solution of (15.31) as combined with (15. 2p is given by 




(5.11) 



where 





For ao > we have n > 0. 



(1) a > 

In this case we can easily check for p > that 




(5.13) 



and, as has been expected, we have 



lim a(p) = ao . 



(5.14) 



(2) a = 

In this case, a = is a fixed point, so that we have a(p) = 0, and consequently 



lim a(p) = . 



(5.15) 



p— »oo 



(3) a < 



15 



The asymptotic limit of a(p), when we start from a negative a, should be equal to either 
or — oo. Since a^o = is repulsive, a^o = — oo is the only choice. In this case we cannot 
use the series expansion in powers of a since \a\ increases indefinitely for increasing p. 



(4) complex a 

Give an imaginary part ie to be added to a 

a ^ a + ie , (5.16) 

then 

l]{l + bg'p)-" + 1^0, (5.17) 



a + te 

so that the formula (15.111) gives 



lima(p) = ao, (e 7^ 0) (5.18) 



p— +00 



for an arbitrary choice of the real part a, provided that the imaginary part is non-zero. As 
mentioned in (1) this result is also valid even for e = for a > 0. The only exceptions are 
the cases a = and a < as mentioned in (2) and (3). 

In Section 4 we have divided the set of real a into three subsets or three classes on the 
real axis, but we can extend this division to the set of complex a as 

=D(i)(_oo)|JZ}W(0)|jD(2)(ao), (5.19) 

where 

D(2) = {a\ aeC}, 
L)(i)(_oo) = {a\ a < 0}, 
D(°)(0) = {a\ a = 0}, 
L)(2)(ap) = {a\ aeC, a<0 excluded}. (5.20) 

The superscripts (2), (1) and (0) denote the dimensionality of these sets or classes, respec- 
tively. 

Next we shall generalize Fig. 1 to the complex plane. The lines of RG flow show their 
behavior qualitatively or topologically, but not quantitatively. 

It is interesting to recognize that the lines of RG flow resemble the lines of force generated 
by a dipole. All lines of flow tend to Oq, except for the one along the negative real axis that 
cannot go up nor down in the complex plane. 
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Figure 3: Lines of RG flow in the complex a plane for increasing p [olq > 0). 



In fl4.29p we found that a = is a discontinuity for the integral of the spectral function 
of the gluon propagator. We can extend this observation to the complex a plane, and for 
this purpose we shall write down the sum rule for negative a and also for a + ie: 

j dm^ p{m^,g^,a,fj,) = 0, (a < 0) (5.21) 

lim / dm^ p(m'^ , g"^ ,a + ie, fi) = — , (a < 0) . (5.22) 
€^0 J ao 

This shows that there is a discontinuity in the integral of the spectral function when a 
approaches the negative real axis. 



QED 

The situation is completely different in QED. First of all, a gauge-invariant two-point 
function 

{FxAx),F^,{y)) (5.23) 

can be expressed in terms of the transverse part of the photon propagator so that the photon 
mass that appears as the pole of this expression is physical. QED is characterized by Ward's 
identity 

Zi = Z2, (5.24) 
which is expressed in terms of anomalous dimensions by 

(3{e) = eMe) , (5.25) 

from which we can derive 

^[e'(p)«(p)]=0. (5.26) 
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As mentioned already, the spectral function p{m?) does not depend on a, nevertheless, we 
have the sum rule: 



The only consistent choice of a^o that makes independent of a is a^o = 0, and = oo 
follows from (15.261) . 

In this case, Ooo is the only asymptotic limit that is attractive, and there is only one 
equivalence class of gauges in QED. Furthermore, 6^0 = 00 signifies that 7y(e) is positive 
definite, and the expression 



is real and decreases with increasing p. Then, the lines of RG flow in QED are shown in Fig. 
4. 



Figure 4: Lines of RG flow in the complex a plane for increasing p (QED). 

Here the lines of RG flow resemble the lines of force generated by a monopole. 
In QED we have only one equivalence class of gauges and the massless photon is physical 
and observable. In this case we have 




(5.27) 




(5.28) 





Z3 = 00 



(5.29) 



and charge confinement is not realized. 
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6 Conclusions 



We may summarize what we have done in this paper in connection with the results obtained 
in a series of papers on this subject [3, EJ [HI fT2] . 

1) Color SU{3) symmetry 

For color confinement we need an exact non-Abelian gauge symmetry. If this symmetry 
is broken, a color singlet state is forced to mix with colored states and consequently color 
cannot be confined. Mathematically, the original form of the condition for color confinement 
was given by the absence of the massless spin zero component in the current SSA'^, but 
breaking of color symmetry would induce the unwanted massless spin zero component in 
that current in the form of the Nambu-Goldstone boson. In Abelian gauge theories such as 
QED the massless spin zero component is always present, thereby preventing the charges 
from being confined. 

2) Condition for Color Confinement 

With the help of renormalization group we can show that a sufficient condition for color 
confinement can be cast in the form of Eq. (14.351) . which formed the basis of our discussion 
in the present paper. 

3) Evaluation of 

Furthermore, with the help of renormalization group we can derive the identity (14.251) 
and can show with the help of asymptotic freedom that for negative choices of a, eventually 
below a certain negative constant, aoo turn out to be —oo, satisfying the condition (14.351) 
for confinement. 

4) ^-matrix in Confined QCD 

By using gauge-independent Green's functions we can express the ^-matrix elements 
for hadronic reactions by applying the reduction formula The unitarity condition for 
the S'-matrix between two hadronic states are expressed by Eq. (I4.36p . and confinement is 
characterized by the statement that the sum over the intermediate states is saturated by 
keeping only hadronic states, excluding colored particles such as quarks and gluons. 

If confinement is realized in one gauge in the above sense, it is clearly realized also 
in all other gauges, since the S'-matrix is gauge- independent. In this sense the concept 
of gauge-independence is a key factor in understanding confinement. The importance of 
gauge-independence or invariance has been stressed by many authors [HI [151 [16] . 
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